Correlation effects on the electronic structure of TiOCl: a NMTO+DMFT study 
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Using the recently developed N-th order muffin-tin orbital-based downfolding technique in combi- 
nation with the Dynamical Mean Field theory, we investigate the electronic properties of the much 
discussed Mott insulator TiOCl in the undimerized phase. Inclusion of correlation effects through 
this approach provides a description of the spectral function into an upper and a lower Hubbard 
band with broad valence states formed out of the orbitally polarized, lower Hubbard band. We find 
that these results are in good agreement with recent photo-emission spectra. 



Introduction - 

A challenging task in the field of strongly electron- 
correlated materials is the description of correlation ef- 
fects, within a possibly controlled approximation, since 
these effects are essential in order to understand the 
behavior of these materials. In the past years there 
has been an enormous effort into improving this de- 
scription by combining ab-initio calculations with many- 
body methods like the LDA+U (Local Density Approx- 
imation plus on-site U)i and the LDA+DMFT (Local 
Density Approximation combined with Dynamical Mean 
Field Theory) 2 . While the original implementation of 
the Dynamical Mean Field theory was based on a single- 
band model, most compounds of interest involve more 
than one correlated orbital. This calls for the need 
of a multi-orbital extension of the LDA+DMFT tech- 
nique. Recently, a new implementation of the multi- 
orbital LDA+DMFT has been proposed which uses the 
localized Wannier functions generated by the N-th or- 
der muffin-tin-orbital (NMTO) method 3 - in order to con- 
struct the LDA Hamiltonian and solves the many-body 
problem by DMFT including the non-diagonal contri- 
bution of the on-site self-energy. Such approach has 
been found to be highly successful in a series of recent 
applications^*^. We will consider this procedure here in 
order to unveil the electronic properties of the layered 
Mott insulator TiOCl. 

TiOCl, which consists of bilayers of Ti 3+ and O 2- par- 
allel to the ab plane, separated by layers of Cl~ ions 
stacked along the c-axis (see Fig.QJ has recently raised a 
lot of discussion due to its puzzling behavior at mod- 
erate to low temperatures. Large phonon anomalies 
have been observed in Raman measurements^ at temper- 
atures around 135K as well as temperature-dependent 
g-factors and line-widths in ESR&. The susceptibility 
measurements^ show a kink at T C 2 = 94 K and an expo- 
nential drop at T c i = 66 K, indicating the opening of a 
spin-gap which is accompanied by a doubling of unit-cell 
along the &-axisi£. While phonons are undoubtedly play- 
ing an important role in the behavior of this system, elec- 
tron correlation is fundamental in order to understand 
these anomalous properties and we will here concentrate 
on the electronic description. 



Recent electronic structure calculations of TiOCl 
within the framework of LDA+U2iii for the crystallo- 
graphic data at T > T C 2 show that in this tempera- 
ture range the ground state of the system is described 
by chains of Ti ions along the 6-direction (see Fig. ^) 
with the 3c? electrons occupying the d xy orbitals. Such 
a study 11 further revealed the importance of the phonon 
degrees of freedom at temperatures T > T C 2 by point- 
ing out that certain A g phonon modes consistent with 
the Pmmn space group of the high temperature struc- 
ture may lead to orbital fluctuations with the ground 
state switching from d xy to d yz /d xz . However, the im- 
plementation of LDA+U calculations 1 — involves the as- 
sumption of a particular spin ordering in the system, 
which is fictitious since there is no true long range order 
in TiOCl. This makes the comparison of the LDA+U 
derived DOS with the experimentally measured photo- 
emission data doubtful. Also, LDA+U being a mean- 
field method, treats the correlation effects beyond LDA 
in a mean-field sense and suppresses all the fluctuation 
effects. For early transition metal oxides like titanates 
which are moderately correlated systems such fluctuation 
effects can turn out to be essentially important. One 
ingenious way of treating the dynamical fluctuation ef- 
fects, is the DMFT which - though it freezes the spatial 
fluctuations- takes fully into account the temporal fluc- 
tuations. In the present communication, we aim to study 
the electronic structure of the high temperature ( T > 
T c2 ) phase of TiOCl by means of the NMTO+DMFT 
method. 

Methodology - 

The octahedral crystal field provided by the O4CI2 oc- 
tahedron, surrounding the Ti site splits the five degener- 
ate d-orbitals into ti g and e g manifolds. The Ti 3+ ion, in 
a 3d 1 configuration, therefore fills 1 /6-th of the t2 g com- 
plex, leaving the e g complex empty. The basic LDA elec- 
tronic structure was reported in RefAii. The essential 
features involve CI and O— p dominated bands extending 
from about -8 to -4 eV, separated by a gap of about 2.5 
eV from the Ti-d complex, with the zero of energy set at 
the Fermi level. The Ti-t2 g bands cross the Fermi level 
with a tiny gap of about 0.2 eV between the occupied t^g 
manifold and unoccupied e g manifold. 
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FIG. 1: Layered structure of TiOCl crystal. Two O atoms 
from the same layer and two O and two CI atoms from neigh- 
boring layers form the O4CI2 octahedra surrounding the Ti 
atom. 

The TiOCl system can therefore be described by a low- 
energy, multi-band Hubbard Hamiltonian, 

H H + — ^ ^ ^mm'^imff^im'-ff (1) 

imra' a 

~\~ ^ ^ ^ (Umm' J rnm 1 )^im(T^iTn'(T 

where n tma = c\ miy c ima and c\ m<y creates a cr-spin elec- 
tron in a localized ti g orbital (m) at site i. H LDA is the 
one-electron Hamiltonian given by LDA. We assume dou- 
ble counting corrections to be orbital-independent within 
the tig manifold, thus resulting into a simple shift of 
the chemical potential. The many-body Hamiltonian 
defined in Eq. 1 depends on the choice of the local- 
ized orbitals. In the present calculation, we employed 
the NMTO^ method for defining these localized orbitals. 
In this method, a basis set of localized orbitals is cons- 
tructed from the exact scattering solutions of a super- 
position of short-ranged, spherically-symmetric potential 
wells (the so-called muffin-tin approximation to the po- 
tential) at a mesh of energies, eo, ei, . . . , e^y. The number 
of energy points, N, defines the order of such muffin-tin 
orbitals, the NMTO's. Each NMTO satisfies a specific 
boundary condition which provides it with an orbital 
character and makes it localized. The NMTO's being 
energy-selective in nature are flexible and may be chosen 
to span selected bands. If these bands are isolated -as 
is in the case of ^-derived bands in TiOCl- the NMTO 
set spans the Hilbert space of the Wannier functions. In 
other words, the orthonormalized NMTO's are the local- 
ized Wannier functions. In our example of TiOCl, the 
low-energy Hamiltonian defined in Eq. 1, involves only 
three correlated, localized ti g Wannier orbitals and no 



other orbitals. However, for such Wannier functions to 
be complete, they must involve contributions from the 
other degrees of freedom, e.g. 0~p and CI— p. This is 
achieved via the downfolding procedure 3 which adds on 
to Ti centered ti g orbitals, tails with O— p and CI— p char- 
acter, thereby defining an effective set of ti g orbitals. It 
is important to note here, that the choice of the orbital 
symmetries d xy , d yz , . . . , d^zi-i depends on the choice of 
the co-ordinate system. Throughout this paper and also 
in previous workSiii, the choice has been made as z — 
a, and x and y axes rotated 45° with respect to b and c. 
However, due to the distorted geometry of the TiOziC^ 
octahedra, the chosen co-ordinate system differs from the 
local co-ordinate system and the various d orbitals de- 
fined following the above mentioned co-ordinate system 
are not the exact eigenstates of the on-site Ti— d Hamil- 
tonian. This introduces mixing between the ti g and e g 
symmetries. The ti g NMTO's therefore, contain also the 
on-site and tail e g character. For a realistic description 
of the LDA contribution and for the proper description 
of the localized orbitals in Eq. 1 it is essential to con- 
sider these issues. To our knowledge, so far no other 
LDA+DMFT is capable of taking these contributions 
into account properly. The plot of such Wannier func- 
tions is shown in Refiii. 

In Fig. [3 we show the off-diagonal part of the 3x3 den- 
sity matrix constructed out of the effective, downfoldcd 
t 2g NMTO's: 

N mm , =J2^ J(E - E n ,k)u™ n 

k.n 

m, m'=l,2,3 refer to t 2g orbitals at the same site and 
u™ n are the appropriate normalized eigenvectors for the 
downfolded t 2g bands in the orthonormalized NMTO ba- 
sis. The summation over k runs over all the k— points in 
the Brillouin zone (BZ) and n runs over all the ti g bands. 
As is seen from the plot, it has non-negligible off-diagonal 
elements within the chosen co-ordinate system. 

The DMFT 2 maps the many-body crystal problem de- 
fined in Eq. 1 onto an effective, self-consistent, multi- 
orbital quantum impurity problem. The corresponding 
local Green's function matrix is calculated via the BZ 
integration, 

GK) = YM"n + H)I - H LDA (k) - EK)]" 1 (2) 

k 

where fj, is the chemical potential defined self- 
consistently through the total number of electrons, Lu n = 
(2n + 1)tt/(3 are the Matsubara frequencies with f3 as 
the inverse temperature ((3 = 1/T). S is the self-energy 
matrix. 

In our DMFT implementation we consider all compo- 
nents of the self-energy matrix S mm ' between different 
tig Wannier functions at a given Ti site, including the 
off-diagonal contribution. The importance of these off- 
diagonal elements has already been inferred from the 
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FIG. 2: Off-diagonal matrix elements of t2 a LDA DOS matrix 
(states/eV) in the Wannier representation. 



plot shown in Fig. 2. This is in contrast to earlier 
LDA+DMFT implementations^ where it was assumed 
that the on-site block of the single-particle Green's func- 
tion is diagonal in the space of the localized orbitals 
defining the many-body Hamiltonian. While such an ap- 
proach is exact for the materials with perfect cubic sym- 
metry, it becomes an approximation for materials having 
distortions. 

The self-consistency condition within DMFT implies 
the local Green's function to be the same as the cor- 
responding solution of the quantum impurity problem 
G (y (t - t') = l/Z/L>[c,ct]e- s ^/c(T)ct(r'), where the 
effective action S e ft is defined in terms of the so-called 
bath Green's function t/~ 1 (w„) = G~ l {w n ) + S CT (o; n ), 
which describes the energy, orbital, spin and tempera- 
ture dependent interaction of a particular site with the 
rest of the medium, 

rP [0 

Seff = - / dr dT'Tr^i^g-^Ty^ir 1 )] 
Jo Jo 

"F ~ ^ ^ U mm' ^ima^Hm' — (j 
imm' tr 

+ ^ ^ ^ {Umm' Jmm' j^ima-IT'im'cr 
im{^m')a 

c(t) = [c.i mo -(r)] is the super- vector of the Grassman 
variables and Z is the partition function. The multi- 
orbital quantum impurity problem is solved by the nu- 
merically exact Quantum Monte Carlo (QMC) scheme. 
Standard parametrization^ has been used for the direct 
and exchange terms of the screened Coulomb interaction, 
Umm' and J mm ' with Umm = U, U mm ' — U - 2J and 
Jmm'(^m) = J- For our calculation we have used U = 4 
eV and J — 0.5 eV, which are reasonable choices for an 
early transition metal like Tiii. The computational ef- 
fort becomes prohibitive rather quickly as one lowers the 
temperature since in order to maintain the accuracy of 
the calculation one needs to increase the imaginary time 
slices as one increases /3. The results reported in the fol- 



lowing are done for {3 = 20 (T = 580 K) with 100 slices 
in imaginary time and 10 6 QMC sweeps. The maximum 
entropy method^ has been used for analytical continua- 
tion of the diagonal part of the Green's function matrix 
to the real energy axis to get the DMFT spectral density. 

In this context, recently a LDA+DMFT calculation 
has been reported for TiOCli&. The calculation scheme 
in Ref»i& is based on a multi-orbital generalization of the 
iterated perturbation theory (IPT) for impurity solver. 
We consider the NMTO+DMFT scheme to be superior 
to the IPT approximation for the present problem. The 
multi-orbital QMC solver gives an accurate solution of 
the correlated orbitally polarized nonmagnetic t2 g prob- 
lem, while the IPT scheme, which is very successful for 
the one-band DMFT problem 2 - becomes quite uncertain 
in the anisotropic multi-band casein. 

LDA +DMFT results and discussions - 

Fig. shows the orbital-resolved and total spec- 
tral functions obtained using the above explained 
LDA+DMFT technique, in comparison with the LDA to- 
tal DOS. Consideration of the correlation effect beyond 
LDA within the framework of DMFT opens up a gap of 
~ 0.3 eV between the occupied and unoccupied spectra, 
signaling the insulating nature of the compound while 
the LDA DOS shows the finite density of states at the 
Fermi energy. We also note the appreciable broaden- 
ing of the bandwidth in LDA+DMFT compared to the 
LDA bandwidth due to the redistribution of the spec- 
tral weight following the opening of the gap. The oc- 
cupied bandwidth obtained in our LDA+DMFT calcu- 
lation is in good agreement with recent photo-emission 
measurements^ which is about 2.5-3 eV. IPT-based 
LDA+DMFT calculations^ gave a bandwidth of 1.5-2 
eV. 

The on-site matrix elements of the tight-binding rep- 
resentation of the downfolded tig Hamiltonian in the 
NMTO Wannier function basis show that the d xy or- 
bital energy is about 0.4 eV lower than the d xz and d yz 
orbital energy, which are degenerate within the chosen 
co-ordinate system. This splitting is much smaller than 
the total tig bandwidth of the LDA DOS which is about 
2 eV and is also smaller than the individual bandwidths 
of d xy and d yz /d xz . As a consequence, in the LDA DOS 
the d xy orbital is occupied by only 0.49 electron with the 
rest occupying the d xz and d yz orbitals. Turning to the 
LDA+DMFT results, the occupation matrix shows that 
the d xy orbital becomes nearly full in contrast to LDA. 
The d xy orbital is found to contain 0.98 electrons. The 
sharp increase in orbital polarization compared to the 
LDA result is in conformity with similar LDA+DMFT 
studies in LaTi03 and YTiOa 4 , but is in contrast to 
IPT-based LDA+DMFT calculations for TiOCl which 
repor1ji& non-negligible inter-orbital mixing, with 70% of 
the electrons residing in the d xy . 

As mentioned already, the gap value obtained in our 
calculation is about 0.3 eV, while the gap extracted from 
the optical conductivity data is about 1 eV in reflectivity 
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FIG. 3: DMFT spectral function at T = 580^ (thick line) and 
the LDA DOS (thin line), [i = 0. The cUy and d V2 /dzg contri- 
butions are shown with dotted and dashed lines respectively. 
The spectral weight of the broad valence states is contributed 
almost entirely by d xy . 

experiments^ and 2 eV in transmittance experiments^, 
the later being claimed as more sensitive for the gap 
determinationSi. As mentioned earlier, our QMC-DMFT 
calculations were carried out at a temperature, T = 580 
K, while the experimental measurements were carried out 
at room temperature ( ps 300K) and below . Recent 
work on shows the successive filling of the Mott- 

Hubbard gap on raising the temperature. Therefore, the 
use of temperatures higher than room temperature may 
be a cause of the underestimation of the gap value ob- 
tained in the theoretical result. Nevertheless, it is hard 



to justify such a large discrepancy as merely a temper- 
ature effect. The IPT-based DMFT calculations which 
claim2& to achieve significantly lower temperatures than 
the QMC-DMFT calculations, predict for TiOCl (with 
U=3.3 eV and J= 1 eV) a gap which is less than 0.5 
e\4i&. The only way to reproduce a gap value as large as 
ps 2 eV, within the single-site DMFT scheme, seems to 
be to use an unreasonably high U value. We therefore 
suspect that a major contribution to this discrepancy is 
caused by the neglect of the component of the self-energy 
between different sites, i.e. by the single-site approxima- 
tion of the DMFT. The dimerization of the crystal and 
formation of spin-singlet bonds at low temperature below 
the ordering temperature can have a precursor at high 
temperatures in form of a short-range ordering effect, re- 
sulting into formation of dynamical Ti-Ti singlet pairs^i 
Consideration of such effects is beyond the scope of the 
single-site DMFT since one needs to consider at least a 
cluster of 2 Ti sites to capture these features. Experi- 
ence with VO2 6 already indicates that such effects can 
be extremely important in widening up the band gap. 
We intend to consider the cluster DMFT for TiOCl in 
our forthcoming communication. We mention here that 
the phonon effects could be important as well, as already 
found in connection with orbital fluctuations, 11 leading 
to additional smearing of the spectral features which is 
probably seen in the experimental PES datai&. 
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